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ABSTRACT 

We demonstrate that the amount of extra mixing required to fit the observed low 
C/N and 12 C/ 13 C ratios in first giant branch (FGB) stars is also sufficient to explain 
the carbon and nitrogen abundances of Galactic asymptotic/second giant branch (AGB) 
stars. We simulate the effect of extra mixing on the FGB by setting the composition of 
the envelope to that observed in low-mass (M < 2Mq) FGB stars, and then evolve the 
models to the tip of the AGB. The inclusion of FGB extra mixing compositional changes 
has a strong effect on the C and N abundance in our AGB models, leading to composi- 
tions consistent with those measured in Galactic carbon-rich stars. The composition of 
the models is also consistent with C and N abundances measured in mainstream silicon 
carbide (SiC) grains. While our models cover the range of C abundances measured in 
carbon stars in the LMC cluster 1846, we cannot simultaneously match the composi- 
tion of the O and C-rich stars at the same time. A second important result is that our 
models only match the oxygen isotopic composition of K and some M, MS giants, and 
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are not able to match the oxygen composition of carbon-rich AGB stars. By increasing 
the abundance of 16 in the intershell (based on observational evidence) it is possible to 
reproduce the observed trend of increasing 16 0/ 18 and 16 0/ 17 ratios with evolution- 
ary phase. We also find that some Li production takes place during the AGB and that 
Li-rich carbon stars (loge(Li) > 1) can be produced. These models show a correlation 
between increasing Li abundances and C. The models cannot explain the composition 
of the most Li-enriched carbon stars, nor can we produce a Li-rich carbon star if we 
assume extra mixing occurs during the FGB owing to 3 He destruction. We tentatively 
conclude that 1) if extra mixing occurs during the AGB it likely only occurs efficiently 
in low metallicity objects, or when the stars are heavily obscured making spectroscopic 
observations difficult, and 2) that the intershell compositions of AGB stars needs further 
investigation. 

Subject headings: Galaxy: abundances - nuclear reactions, nucleosynthesis, abundances 
- stars: Abundances, stars: AGB and Post-AGB 



Introduction 



Standard stellar modela^ predict that there is only one mixing event that changes the surface 
composition between the main sequence and the tip of the first giant branch (FGB) - the first 
dredge- up (FDU). The models indicate that the FDU mixes the products of partial hydrogen 
burning into the envelope, with reductions in the 12 C/ 13 C ratio to ~ 20, and increases in 14 N such 
that 12 C/ 14 N ~ 1.5. The 16 Q/ 17 Q ratio is reduced, depending on the mass of the star, whereas 
16 0/ 18 marginally increases (e.g., Boothroyd Sackmannlll999 ). Lithium, if not destroyed during 
the pre-main sequence, should be greatly reduced by the FDU. In contrast the envelope abundance 
of 3 He, produced by the pp-chains during the main sequence, is significantly increased by the FDU. 
Observations, however, do not agree with these predictions. Low-mass giant sta r s have been found 
with 12 C/ 13 C ratios of ~ 10 and C/N ~ 1.0 (jGilrov! Il989l : Icharbonnell Il994l : Icharbonnel et al 



19981 ). In globular clusters, the anticorrelation of C abundance with luminosity is also indicative 
of deep mixing occurring on the FG B (e.g., in M3: Smith 20 021). Furthermor e, chemical evolution 
models using standard 3 He yields ( Dearborn et al. 19961 : Romano et al. 2003 ) predict present-day 
abundances that are orders of magnitude greater than the values measured in meteorites and the 
local interstellar medium. These observations are used as evidence that some form of chemical 
transport (known variously as extra-mixing, deep mixing, or cool bottom processing) is acting in 
low-mass red giant envelopes. 



1 Standard in this context refers to models without non-convective extra mixing such as thermohaline mixing or 
rotation. This is distinguished from convective overshoot, which is also non-standard in stellar codes but is considered 
a natural extension of convective borders. 
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In recent years there has been m uch effort devoted to fin ding the physic al mechanism behinc 



the e xtra mixing. Rotational mixing (Charbo nnel et al.lll998f). gravity waves (Dcnisse nkov Tout 



2009: IStancliffdl2010h . and magnetic fields ([Nordhaus et alJl2008l : lBusso et al.ll2007l : IPalmerini et al 



20001 ) thermohaline mixing (|Eggleton et alJbood. bood: Icharbonnel & Zahn)[2 007: 



Stancliffe et al 



20091 ) have been proposed as potential candidates. Parametric studies of extra mixing have also 
shown great promise in explaining the C, N, and Q compos i tion of first and asy mptotic giant branch 



stars ([Smith k, Tout! LL992J; iBoothroyd fc Sackmann 



1999 



Nollett et al.ll2003l ). although there are 



no compelling physical reasons for the choice of parameters that govern the mixing and burning. 

The evidence that some form of non-convective mixing mechanism occurs during the AGB 
phase is more circumstantial. The argument put forward is similar to that for FGB giant stars. 
That is, the C and N abundances predicted by models do not m atch those observed in AGB stars 



(e.g.. lLambert et al.lll986l : lAbia &: Isern 



1997 



Milam et al 



2009) . In particular the 12 C/"C ratios 



Forestini &; Charbonnell 1997 : Karakas 



are l ower than predicted by standard AGB models (e.g. 
2010) . Also, the C and N composition of pre-solar mainstream silicon carbide (SiC) grains, which 
are assumed to form in the extended envelopes of carbon-rich (C/O > 1) AGB stars, show a well 
defined distribution where 40 < 12 C/ 13 C < 100, and have nitrogen isotope ratios enhanced relative 
to the terrestrial value (e.g.. IZinnerlll998l ). In comparison, standard AGB models predict that by 
the time C/O > 1 the 12 C/ 13 C > 80, which is already close to the upper limit observed in AGB 
stars or measured in SiC grains. A few studies have found that the only way to match observational 
data of carbon-rich AG B stars is to artificial l y lower the 12 C/ 13 C a t the tip of the FGB to values 



observed in FGB st ars (jKahane et al.1 12000 : iLebzelter et al.l l2008l ). For example, by doing this 



Kahane et al.1 (j2000i ) found the abundance ratios in their 2M & model of CW Leo were consistent 
with the observations. However, the authors did not go on to test this idea against data for a 
larger sample of Galactic carbon stars. Lithium is also enhanced in some carbon-rich AGB stars. 
The statistics indicate that super Li-rich stars with Li abu ndances as high a s loge(Li) ~ are 
rare, accounting for only a few percent of all carbon stars (jAbia et al.lll993l ). However a larger 
fraction (~ 10%) are considered Li- rich, with Li abu ndances in the range 1-2 (although some 
J-stars are included in this fraction, lAbia et al.l Il993l ). Standard low-mass AGB models do not 
predict this variation in Li, thus extra mixing may be needed to explain the high Li abundance s 
in some AGB stars (jAbia Isernl Il997l ; luttenthaler et al.l bood : luttenthaler fc Lebzelterl bold ). 
Intermediat e-mass AGB stars can produce substantial quantities of Li via hot bottom burning 
(HBB; e.g.. IBoothroyd et al.lll993l ). While most C-rich stars will have a mass below the HBB limit 
(M < 4Mq), a few super-Li rich stars may have evolved from high mass AGB stars. These points 
are discussed further in 

If we assume that thermohaline mixing (or some other mechanism) operates to change the 
envelope composition during the ascent of the FGB, then a natural questions arises: Do we need 



2 We use the standard notation loge(X) = log 10 (Nx/Nn) + 12.0, where iVx is the number abundance of species X 
and Nn is the number abundance of hydrogen. 
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further extra mixing on the AGB to be able to match the observed data? It is this question that 
we seek to address here. This question is particularly timely, given the recent ef forts to find the 
physical mechanism driving the extra mixing on the FGB (e.g.. lEggleton et al.ll2008l ). In this paper 
we argue that there is no strong need for extra mixing on the AGB, at least for Galactic carbon 
stars of near solar metallicity. For the moderately metal-poor AGB s tars in the LMC and SMC 
the situation is less clear (see also the discussion in lLederer et alJbood ). 

In the current study we evolve a series of low-mass AGB models, some of which we assume 
have experienced extra mixing on the FGB. We then compare our models with observations of C/O 
ratios and carbon isotopic data in a sample of Galactic carbon stars and mainstream SiC grains, 
and the compositions of stars in the Magellanic Clusters NGC 1846 and NGC 1978. We include a 
discussion of lithium and the oxygen isotope ratios in AGB stars, as these species may elucidate 
the occurrence and nature of possible extra mixing in AGB stars. 



Stellar Models 



We evolve 1.25M and I.8M0, Z = 0.01 stellar models from the zero age main sequence 
(ZAMS) to the tip of the AGB. These masses were chosen because the 1.25M is likely to be close 
to the lower mass limit for carbon star format i on, whereas the 1 .8Mp) is a typical carbon star mass in 



both the Galaxy (jWallerstein Knapp 



1998 



Abia et al 



2001 



and in the LMC (Lebzelter fc Wood 



20071 ). In regard to the metallicities of Galactic carbon stars, I Lambert et al.l (j 19861 ) noted that 
—0.3 < I'M/HI < 0.0, where M is the total metal abundance. These results were confirmed by 
Abia et al.l (|2002l ) who confirm that the metallicities of N-type carbon stars are of nearly solar 
metallicity. These results suggest our choice of Z = 0.01 (or [Fe/H] = —0.15) is consistent with 
the metallicities of Galactic carbon stars. We calculate the structure first and perform detailed 
nucleosynthesis calculations afterward, using a post-processing algorithm. The details of this pro- 
cedure and the cod es used to compute the models have been previously described in detail (e.g., 
Karakas et al J 12009 ). In the nucleosynthesis step we use 77 spec i es, an d a scaled-solar initial com- 
position where we take the solar abundances from lAsplund et al.l (|2005l ) and assume Z S oiar = 0.015. 
Mass loss on the FGB is included usin g the R eimer's formula with the parameter 77 = 0.4 and on the 
AGB we use the IVassiliadis fc Wood! ([ 19931 . hereafter VW93) mass-loss prescription. The VW93 
mass-loss prescription keeps the mass-loss rate relatively low until the beginning of the superwind 
phase, at which point the mass-loss rate dramatically increases to a few xlO" 5 M year . This 
has the effect of keeping the envelope mass relatively high, compared to using the Reimer's mass- 
l oss r ate on the AGB, and may affect the depth of third dredge-up (TDU). IStancliffe fc Jeffery 
(|2007l ) discuss the effect of varying the mass loss on the evolution and nucleosynthesis of a 1.5M , 
Z = 0.008 model. They find that the total amount of material that gets dredged- up to be similar 
for all three mass-loss prescriptions (VW93, Reimer's with rj = 1, and [Blocker! 1 1 995l . with 7/ = 0.02). 
The surface C/O ratio, however, tended to be higher at a given time in the Reimer's model than in 
the other two cases. The reason given for this is because models using Reimer's law typically had 
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a small envelope mass at a given thermal pulse number, It should be noted that both Reimer's and 
the Blocker mass loss critically depend on the choice of the free parameter rj, the values of which 
are uncertain. 

The main changes to the stellar structure code that have been imp lemented for this study ar e 
the inclusion of the C and N-rich low temperature opacity tables from lLederer &: Aringerl ([2009), 
and the inclusion of convective overshoot to induce the TDU in low-mass models. The TDU is the 
inward movement of the convective eny elope into regions processed by partial He-burning during 
a thermal instability (see iHerwia 120051 . for a recent review) . Some form of overshoot is required 
because low-mass models computed prev iously with the sam e codes show little or no TDU for 
masses < 2M & at near solar metallicities (jKarakas et al.ll2002l ). Other stell ar evolution codes also 



show no TDU at these masses without some form of convective overshoot ( Mowlavi 



1999 



Herwig 



2000). We include overshoot following a thermal pulse by extending the position of the base of 
the convective envelope by N ov pressure-scale heights. To obtain a carbon-rich envelope we require 
0.5 < N ov i$ 5, depending on the mass of the model. In the 1.25M model we set N ov = 3,4, 
and in the 1.8M Q model we set N ov = 1,2,3. Experimenting with the efficiency of overshoot in 
AGB models is reasonable, given our lack of understanding of the physics of convection in stars. 
The location of the inner edge of the convective envelope has been sho wn to depend critically o n 
the numerical details of a particular stellar evolution code (see Fig. 1 in iFrost h Lattanzid Il996l ) , 
and d ifferent stellar evolution codes predic t dissimilar results fo r the same masses and metallicities 



(e.g., IStraniero et al.lll997l : iMowlavi 



1999 



Karakas et al 



2002 



Stancliffe et al.1 12004). Efforts to 



match the carbon-star luminosity function in the Magellanic Clouds ty pically assu me more efficient 



TDU than found in low- mass stellar models (e.g., iMarigo et al.lll999l . but see also IStancliffe et al 
20051 .) Furthermore, it is still unknown if the efficiency of TDU is stochastic in nature, or there is 
one depth, for a given total stellar mass, H-exhausted core mass, and composition. 

In Table [T] we show the initial, post-FDU, tip of the FGB for the extra- mixing cas^ 3 ], and tip of 
the AGB ratios for C/O, 12 C/ 13 C, 12 C/ 14 N, 14 N/ 15 N, 16 0/ 17 0, and 16 0/ 18 from the stellar mod- 
els. The post-FDU values of 12 C/ 13 C and C/N for the 1.25M© model are 27 and 1.48, respectively. 
The post-FDU values for the 1.8 Mp) model are 22 and 1.0. These numbers are consistent wi th 
other standard model predictions (IE1 Eidlll994l : ICharbonnellll994l : iBoothroyd &; Sackmannlll999l ). 



To simulate the effect of extra mixing on the FGB, we take the envelope compositi on at the tip 
of th e giant branch and alter it such that the 12 C/ 13 C ratio equa ls 10 (from Fig. 6 of ICharbonnel 



19941 ) , and the C /N ratio ~ 1 (Table 1 in iGilrov &; Brownl Il99ll ) . This is done by decreasing the 
12 C abundance and increasing the abundances of 13 C and N, as expected from CN cycling. The 
composition of the oxygen isotopes are not altered because the observed data for FGB stars are 
consistent with model predictions (jDearbornl 119921 : IBoothroyd <fe Sackmannl Il999l ) . By using the 
lower limit for 12 C/ 13 C from the observations (which show values between 10 and 18) we expect to 
maximize the effect on the final AGB composition. We would expect intermediate values to give 



3 for the models without extra mixing, the abundances at the tip of the FGB are the same as after the FDU. 
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results that fall between our extreme case and the no mixing model. For the 1.8M & case the C/N 
ratio is already 1.0 after FDU. After "extra-mixing", a reduction in 12 C and increases in 13 C and 
14 N result in C/N « 0.9 (see Table [1]). This value is still within the range of observed C / N ratios 
observed in FGB stars, that show C/N ratios between ~ 0.8 to 1.2 (jGilrov Brownlll99ll ). 



Along with simulating the effect of extra mixing on the composition of the envelope, we perform 
two further tests. The first is to increase the 16 abundance at the ZAMS so that the initial C/O 
ratio = 0.35 (in comparison to other models which have a solar C/O = 0.5). This is equivalent 
to an a-enhanced composition, where [O/Fe] ~ +0.20. We do this because F and G stars in our 



Galaxy show moderate enhancements of a elements (e.g.. lReddv et al.ll2003l ). and perhaps some of 
the Galactic AGB stars have similar moderate enhancements. 

The second test is to increase the intershell abundance of 12 C and 16 to values observed in 
the central stars of planetary nebulae. Standard model predictions provide values of about 25% 
(by mass) of 12 C in the intershell, and only < 2% of 16 0. Observations of post-AG B central stars 



sugge st the intershell abundances are between 2-20% for 16 and 15-60% for 12 C (jWerner et al 
20091 ). 



3. Model predictions versus real stars 



In Fig. [I] we show the distribution of C/ C ratios from lLambert et al.l (j 19861 ). The sample 



average C/ C ratio for the C-rich stars is 58, excluding the J-stars, or 50 with the J-stars 
included. The standard deviation on these numbers are 25 and 17, with and without the J-stars, 
respectively. We choose to use this data-set for comparison because these authors use infra-red 
observations to derive the 12 C/ 13 C ratios and are likely more reliable for the very red, cool carbon- 
rich AGB stars under consideration here. We exclude the 12 C/ 13 C ratios for the J-type carbon 
stars, which have carbon isotope ratios close to the equilibrium value of ~ 4. The evolutionary 
origin of these 13 C-rich stars has not been determined. Many J-stars do not show enrichments in 
g-proc ess elements, in contrast to N-type stars and it is questionable if they are really TP-AGB 
stars ( Wallerstein &: Knapplll998l ). The average value of the 12 C/ 13 C ratio for main stream SiC 
grains is 57 with a standard deviation a = 16 (using data from iGyngard et al.ll2006l ). consistent 
with the average value derived for the C-rich stars in Fig. [TJ 

Included in Fig. Q]are model predictions from the 1.25Mq (top panel), and the 1.8M (lower 
panel) AGB models. In the top panel, the dashed vertical line shows the 12 C/ 13 C ratio at the tip 
of the AGB for the N ov = 3 model with no extra mixing on the FGB, whereas the dot-dashed 
vertical line shows the result when the 12 C/ 13 C ratio is reduced at the tip of the FGB in line with 
observations. This clearly shows the strong effect on the carbon isotope ratio in AGB stars when 
taking FGB mixing into account. In this example the 12 C/ 13 C ratio is reduced from ~ 76 to 33. 
The dotted vertical line shows the 12 C/ 13 C ratio for the N ov = 4 case with extra mixing on the 
FGB. This shows the effect of deeper TDU increasing the isotopic ratio due to the extra 12 C that 
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is dredged up. Likewise for the bottom panel the dashed vertical line shows the 12 C/ 13 C ratio at 
the tip of the AGB for the N ov = 1 model with no extra mixing, whereas the results including 
extra mixing are shown by the vertical dot-dashed line. Again, we can see the significant effect 
of including the FGB abundance ratios. The dotted vertical line shows the 12 C/ 13 C ratio for the 
N ov = 2 case with extra mixing (without FGB extra mixing this model has a 12 C/ 13 C ratio > 100, 
off the graph) . In summary, from Fig. Q] we see that the models considered here span the observed 
range of 12 C/ 13 C ratios, with th e exception of one ob ject which has a ratio of 19 ± 6. Two objects 
have 12 C/ 13 C ratios ~ 30 ± 5 ( Lambert et al.l ll986) but the uncertainties (±5) allow these two 
stars to be spanned by the model predictions. Furthermore, it is possible that not all C-stars will 
be at the tip of the AGB. Some of the carbon stars may have a value between the FGB tip (after 
modification by extra mixing) and the AGB tip value. 

In Fig. [2] we illustrate the temporal evolution of the C/O ratio versus the 12 C/ 13 C ratio for 
the 1.8M model with N ov = 3 (open symbols and lines, refer to the caption for details). Note 
that each symbol represents the model abundances after a TDU episode. This model was chosen 
because it reaches a high final C/O ratio of 3.3 (models with different values of N ov follow similar 
tracks, but finish at lower C/O ratios). It will become apparent that once we start enhancing the 
O abundance (either at the beginning of the evolution, or in the intershell) the final C/O ratios will 
be reduced and for this comparison we require the final C/O ratio > 1 . In Fig. [2] we show results 
from models with scaled-solar and a-enhanced compositions, with and without extra mixing on the 
FGB. The labels on each line show the C/O and 12 C/ 13 C ratios at the tip of the FGB. 

Note that the models with an a-enhanced composition (i.e., lower initial C/O ratio) shift the 
predictions to the left, and result in a steeper slope and low er final C/O rat i o. A lso included in 
Fig. [2] are observational data for Galactic carbon stars from lLambert et al.l (1198m. and s t ars in 
the Magellanic Cloud clu sters NGC 1846 and NGC 1978 using data from iLebzelter et aD (|2008h 
and iLederer et al.l (|2009l ) . respectively. These clusters have a lower metallicity ([Fe/H] ~ —0.5) 
than considered her e ([Fe/H] ~ —0.15), and the AGB mass of NGC 1978 is smaller, at 1.5Mq 
(jLederer et al.1 l2009l ) . It is worth noting that the abundance tracks of the 1.25M0 model (not 
shown) follow almost the same path as the 1.8M & models, the only noticeable difference being a 
final lower C/O ratio as a consequence of fewer thermal pulses and TDU episodes. Even so, it 
is clear from Fig. [2] that the various stellar models cover most of the range of observational data 
points. This agreement is especially encouraging considering that we are using only one mass and 
one metallicity to cover the range of observed carbon-rich AGB stars, where there is likely a spread 
in initial mass, composition, and evolutionary state. 



The hypothesis put forward by ILebzelter et al.l (|2008l ) to explain the composition of the NGC 
1846 stars is as follows. The O-rich stars were matched by assuming extra mixing on the FGB 
plus an a-enhanced composition ([O/Fe] = +0.2), such that the C/O ratio = 0.2 at the tip of the 
FGB, consistent with the observations which yield C/O = 0.2 ± 0.05. To explain the composition 
of the two carbon-rich NGC 1846 stars near the bottom right of Fig. [2] required extra mixing on 
the FGB and AGB, but mixing only occurs on the AGB after the stars had become carbon rich. 
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It is not easy to understand why mixing would only start once C/O > 1 in the envelope, given 
that there are several episodes of TDU before this. While the outer envelope (e.g., the opacities) 
strongly depend on the composition, the structure of the star near the inner region of the envelope 
and the H-shell do not. There are no noticeable changes to the structure of the burning shells, nor 
in the time between pulses as the star evolves from C/O < 1 to carbon rich. Fig. [2] shows that 
our scaled-solar models with FGB extra mixing alone can explain the composition of these two 
C-rich stars. However, the models that fit the t wo stars have a higher (0.32) C/O ratio at the tip 
of the FGB than observed in the O-rich stars by iLebzelter et al.1 (|2008l ) . Hence our models cannot 
simultaneously fit the composition of the O and C-rich stars at the same time. One test of the 
AGB extra mixing scenario would be the determination of the Li abundances in the two stars. A 
high Li abundance may imply that these stars have experienced extra mixing on the AGB whereas 
a low Li abundance may rule out the need for extra mixing. 

Lastly in this section we discuss the results for lithium from our stellar models. We will use the 
1.8Mq model with A ov = 3 as an example. While some C and Li-rich stars with C/O ratios ~ 1 
may have evolved from high-mass AGB stars undergoing HBB, most C-rich objects evolved from 
lower mass stars. The first dredge- up reduces the Li envelope abundance from loge(Li) from 3.14 
at the ZAMS to 1.28. It is likely that our models do not suffer enough Li depletion from the first 



dredg e-up, where observations of G-K giants suggest most stars have loge(Li) < 1.0 (jBrown et al 



1989). During the TP- AGB the 7 Li abundance in our models increases as a consequence of the 
TDU, and reaches a maximum of loge(Li) = 1.80. The increase of Li mostly takes place before 
efficient TDU begins to penetrate into the He-intershell (i.e., during the first few TPs which are 
weak), but after this we find a correlation between the Li abundance and carbon (and presumably, 
g-process elements). This could explain the correlation f ound between Li and Tc abundances by 
Vanture et al. (j2007l ) and Uttenthaler fc Lebzelter ( 2010 ). We note that the final Li abund ance 
of this model would be considered Li-rich according to the definition of lAbia et al.l ()1993l ) i.e., 
loge(Li) > 1.0. If our model had suffered extra depletion on the FGB, then it is possible that the 
Li abundance would not have reached log e(Li) > 1.0 during the TP- AGB. We test this by reducing 
the envelope 7 Li abundance at the tip of the FGB to loge(Li) = 0.0 (close to the center of the Li 
abundance distribution for G-K giants, Fig. 3 from lAbia et al.lll993l ). We also reduce the envelope 
3 He abundance by a factor of 10, because this isotope would also be destroyed by whatever process 
destroys the lithium. While we find a similar trend as before, with Li increasing prior to the start 
of deep TDU and then slowly thereafter, the final value is 0.76, less than the minimum for Li-rich 
carbon stars. 



3.1. The oxygen isotope ratios 

In this section we compare our models to oxygen isotopic measurements of FGB and AGB 
stars. In particular, we are interested to see if standard AGB models can match the observed 
composition, and if not, do we need to resort to extra mixing on the AGB? We focus on oxygen 
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isotopic data from the stellar spectra of evolved stars and do not compare to the oxygen composition 
of presolar oxide grains. We do this because the oxygen isotope data for stars spans a smaller range 
than measured in grains, with few stars found in the region occupie d by Group 2 grain s, and no 



stars found in the region occupied by Groups 3 and 4 (Fig. 9 from iNittler et al.lll997l ). If these 
grains formed in the outflows of AGB stars then there should be stellar counterparts. Or perhaps 
the grains formed in circumstellar envelopes around AGB stars much more evolved (and heavily 
obscured) than the MS, S and C-type giants observed by Harris and collaborators. If so, it is 
curious that the C isotope data measured in SiC grains span the same range as measured in AGB 
stars while the O isotopic data in oxide grains do not. 

In Fig. [3]we show the 16 0/ 18 ratio versus the 16 0/ 17 ratio from a sample of G and K giants, 
barium stars, and carbon-rich AGB stars. References are provided in the figure caption. Error bars 
are not included but uncertainties are large; we refer to Harris et al. (jl987l ) for a detailed discussion. 



Included in the figure are our predictions from the 1.25M and 1.8M@ scaled-solar models. We 
note that we are using models without extra mixing on the FGB, but the results would not change 
because we assume that the oxygen isotope ratios are not altered by this mixing. We also indicate 
in the figure the solar oxygen isotope ratios, which provide the starting position for our models 
(the dotted lines). In Tableland Fig. [3j it can be seen that the FDU reduces the 16 0/ 17 ratio 
to ~ 330 in the 1.8M & model, and 1800 in the 1.25M case, as a consequence of CNO cycling that 
synthesizes some 17 0. In contrast, the 16 0/ 18 ratio increases slightly as a consequence of 18 
destr uction. These numbers a re reasonably consistent with the O isotope ratios measured in FGB 



stars (jHarris fc Lambertlll984l ) 



It is also evident in Table [JJ and Fig. [3] that the predicted oxygen isotope ratios at the tip of the 
AGB do not match the isotopes measured in C-rich AGB stars. In particular, the oxygen isotope 
ratios derived from observations of AGB stars are notably larger than the model predictions. Is 
this caused by AGB extra mixing? This is unlikely, as the changes caused by extra mixing would 
shift the 16 0/ 17 and 16 0/ 18 ratios in a similar direction as the FDU. That is, the abundance 
of 17 would be enhanced, 18 destroyed, and 16 would remain largely unchanged. This would 
result in increases in 16 Q/ 18 Q and decreas es in 16 0/ 17 0. Instead, both ratios appear to increase 



in the stars with evolution. iHerwigl (|2000l ) also noted that the oxygen isotope ratios measured in 
AGB stars are correlated with each other, and seem to be correlated with the observed carbon and 
s-process element abundances. Hence it is likely that the increase of the oxygen isotope ratios with 
evolutionary state can be attributed to the TDU. 

We now explore another possibility to explain the discrepancy. The operation of the TDU 
mixes material from the He-burning shell, that has been exposed to partial He-burning, to the 
surface. The composition of the intershell is mostly 4 He and 12 C (with typical compositions 75% 
and 25%, respectively). Partial He-burning also synthesizes small quantities of 16 and 22 Ne, at 
the level of a few percent each at most. Could it be that the amount of 16 O in the intershell is 
higher than predicted by standard models? 
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Support for this idea comes from both observational and theoretical fronts. The composition 
of hot post-AGB stars show C abundances between 15-60%, O abundances between 2-20%, an d 
Ne abundances of about 2% (be Marco & Barlowl boOll : [Wesson et al.ll2008l : Iwerner et al.1 l2009h . 
The theoretical models of iHerwia (|2000l ) with exponentially decaying overshoot at all convective 
borders find O abundances in the intershell of ~ 20%. The overshoot penetrates into the top of 
the C-0 core, increasing the C and O content of the intershell during a convective thermal pulse. 
Furthermore, key He-burning reactions rates (e.g., 12 C(a,7) 16 0) are uncertain and could allow for 
higher 16 intershell abundances. 



The next question is how do these intershell compositions effect the envelope composition? 
To test this idea, we synthetically compute the envelope abundances of the 1.25Mq model with 
N OY = 4, and 1.8Mq model with N ov = 3. All tests assume a scaled-solar initial composition and 
that extra mixing has taken place on the FGB. From the detailed stellar structure computations, 
we have the total mass, H-exhausted core mass, and mass dredged into the envelope as a function 
of thermal pulse number. Assuming a given set of intershell abundances we then calculate the 
envelope abundance as a function of pulse number. In Table [2] we show the 12 C and 16 intershell 
abundances used (in mass fractions), and the resultant surface abundances at the tip of the AGB 
for the C/O, 12 C/ 13 C, 16 0/ 17 0, and 16 0/ 18 ratios. The surface nitrogen isotope ratio is assumed 
not to change as a consequence of increased oxygen in the intershell. These results are illustrated 
in Fig. [3] for the oxygen isotope ratios, and in Fig. U] for the C/O and 12 C/ 13 C ratios. In Fig. [3] 
the open squares show the result of using an intershell composition of 20% for O in the 1.8Mq 
model. Note that the oxygen isotope ratios do not significantly shift from the post-FDU values for 
the 1.8M model assuming a standard intershell composition. The open circles show the evolution 
of the O isotope ratios in the 1.25M & case with a standard intershell composition (26% C and 
0.08% O), whereas the open triangles show the result of using 20% O. The two models with O-rich 
intershells show a similar trend of increasing oxygen isotope ratios with evolution as seen in the 
observations. 

The comparison is not perfect, however, but the tests performed here are simple and rely 
on many assumptions. First, we have assumed constant intershell abundances, whereas detailed 
models show that the He-intershell abundances of 12 C and 16 O depend on the total mass, core mass, 
initial composition, and time during the convective thermal pulse. Second, we have only enhanced 
the abundance of 12 C and 16 0, and ignored changes to the minor 17 and 18 isotopes. This is 
a valid assumption given that the C-0 core below the He shell contains only trace quantities of 
these isotopes. Also, we have only considered models of two masses, whereas real carbon-rich AGB 
stars in the Galaxy likely evolve from initial masses between about 1.2M Q to ~ 3Mq. Certainly 
models of stars of mass between 1.25M0 and 1.8M & may help fill in the gap between the model 
predictions in Fig. [3j We also do not vary the rates of the uncertain l7 0(p, 7) 18 F and l7 0(p, a) 14 N 
reactions, which play an important role in setting the minimum 16 O/ 1 O value reached after the 



FDU. The new evaluations of the 17 +p rates by IChafa et al.1 ( 20071 ) have reduced uncertainties 
compared to previous estimates, but these still range between ~ 20% to a factor of two, depending 
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on temperature. 



In Fig. U]we show the results of an increasing C and O intershell composition on the evolution 
of the C/O and 12 C/ 13 C ratios from the 1.8M model with N ov = 3. The solid line shows model 
predictions assuming a standard intershell composition (24% C and 0.6% O) and extra mixing on 
the FGB. The dashed and dot-dashed lines show the result of increasing the 16 O content of the 
intershell to 10% (dashed) and 15% (dot-dashed), which has the effect of shifting the final C/O 
ratio to lower values. We note that the extra 16 from the intershell helps to keep the final C/O 
ratio within the observed range. The dott ed line illustrates results that use the composition of PG 
1159-035, which has 48% C and 17% O (jjahn et al.l 120071 : IWesson et alJl^OQsl ). in the intershell. 
That the model predictions in Fig. 2] fit the data for Galactic carbon stars remarkably well may 
indicate that current AGB models predict too little 16 in the intershell. This uncertainty warrants 
further investigation. 



4. Discussion and concluding remarks 



We now return to the original question posed in the Introduction: Do we need extra mixing 
on the AGB? If we simply seek to explain the C/O and 12 C/ 13 C ratios in Galactic AGB stars 
then we have found the answer to be no. The inclusion of extra mixing on the FGB followed by 
a standard thermally-pulsing AGB evolution can produce C and N compositions consistent with 
values observed in AGB stars and measured in pre-solar SiC grains. 

But what about Li in carbon-rich Galactic AGB stars? The existence of stars that are both 
C- and Li-rich has led some autho rs to conclude that extra mixing occurs in some fraction of 
AGB stars (e.g., lAbia &: Isernlll997l ). The frequency of carbon stars rich in Li is < 10%, implying 
that the majority of carbon stars are Li normal or Li poor (loge(Li) < 1). Furthermore, there is 
evidence that the Li production is tied to thermal pulses (and possibly the TDU) in AGB stars, as 
indicated b y the correlation between AGB stars rich in Li an d enriched in the unstable s-process 
element Tc (jVanture et al.ll2007l ; lUttenthaler h Lebzelterll2010l ) . We found that lithium is produced 
during the TP- AGB phase of the 1.8M model with deep TDU, and we find a correlation between 
increasing Li and carbon in the envelope (and presumably also s-process ele ments) . That our models 
can re produce the correlation between He-shell products and Li found by lUttenthaler &: Lebzelter 
(|2010l ) is compelling, but further study into Li nucleosynthesis is required. Our model can only 
become a Li-rich carbon star when we assume that no extra mixing has taken place during the FGB. 
This model then has a higher 12 C/ 13 C and C/O ratio than observed. If we assume extra mixing has 
taken place on the FGB, the star still produces Li but is no longer considered "Li-rich", hence extra 
mixing would be required during the AGB to provide the Li. However we restate that only ~ 10% of 
carbon stars are observed to be Li-rich, which indicates that efficient extra mixing during the AGB 
is not required in the majority of these objects. We note that this is consistent with our findings 
for C and N. The spread in Li observed in carbon-rich AGB stars may be explained by the fact 
we have used one stellar mass and metallicity for our tests whereas Galactic carbon stars will have 
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a range of stellar masses, s ome of which are pre dicted not to experience any extra mixing during 
the FGB (i.e., M > 2M & lEggleton et al.l l2008h . Also, it is possible that stars experience some 
variation in the mixing efficiency, with the Li abundance able t o increase or decrease depending on 
the details of the mixing and burning (jWasserburg et al.lll995l ). 



The low C/O ratios present in C-rich AGB stars as compared to theoretical models (see e.g., 
Fig. [2]) have been previously explained by allowing for extra mixing on the AGB or by the use of C- 
rich low-temperature opacities. Extra mixi ng would reduce the C abund a nce b ut keep the elemental 
O abundance approximately constant (e.g- jNollett et al.ll2003l ). iMarigol ([20021 ) showed that the use 
of C-rich low-temperature opacities leads to an increase of the mass-loss rate and therefore limits 
the final C/O of the model star. I n our detailed models we use the latest low-temperature C-rich 
opacities from lLederer Aringerl (120091 ) and we still find higher final C/O ratios than spanned by 
the observations, although these models depend on the assumed efficiency of convective overshoot. 
Even so, our models suggest that there are only two ways to reduce the C/O ratio in AGB models: 
Variations in the intershell abundances or extra mixing on the AGB. In Fig. H]we show that models 
with high O intershell abundances are able to produce C/O ratios within the range measured in 
Galactic carbon stars. With a reduced TDU, these same models may also be able to account for 
AGB stars rich in s-process elements but with spectral types M and MS, which have C/O < 1. We 
have shown extra m i xing n eed not be the only explanation for these observations. This was also 
discussed by iHerwia (|2000l ) although his models were more massive (3,4M Q ) than considered here 
and likely more massive than most carbon stars. 

The case for extra mixing on the AGB is much stronger at low metallicity. The LMC clusters 
NGC 1846 and NGC 1978 both have a m etallicity lower than sol ar, although NGC 1846 is slightly 
more met al poor with [Fe/H] = -0.49 (|Grocholski et al.l bood ). whereas NGC 1978 has [Fe/H] 
= —0.37 (jMucciarelli et al.l 120071 ) . To explain the composition of the O and C-rich stars in NGC 
1846, extra mixing was proposed to occur on the FGB and AGB, but only on the AGB once the 
stars had become carbon rich. In contrast, extra mixing on the AGB does not seem to be required in 



the L MC cluster NGC 1978, where the C-rich stars have high 12 C/ 13 C ratios (> 150) (jLederer et al 



2009). Why the evolution of the C-rich stars in the cluster NGC 1978 should be so different to 
those in the cluster NGC 1846 is puzzling, given the small difference in metallicity. From the trend 
of the models in Fig. UJ we speculate that an intershell enhanced in both 12 C and 16 may be 
able to explain the abundances of the LMC cluster NGC 1978. Stellar models of the appropriate 
mass and metallicity as the AGB stars in that cluster are needed to explore this possibility in more 
detail. 

At lower metallicities ([Fe/H] 5~ —2) very low 12 C/ 13 C r atios are measured in carbon-enhanced 



enrichments (Aoki et al. 2008 



metal-poor ( CEMP) stars (e . g. . iBeers &: Ch ristliebj 20051 : ISneden et al.1 120081 ). as are high Li 



Thompson et al.l 120081 ) . The origin of the CEMP stars with en- 



richments of s-process elements has been suggested to be mass transfer from a previous AGB 
companion. This is supported by observa tions that suggest th at CEMP stars with s-process ele- 
ments are all members of binary systems (jLucatello et al.ll2005l ). Extra mixing on the FGB of the 
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donor AGB star would not be enough to explain the low carbon iso topic compositions measured in 
the CEMP stars (jSivarani et al J 120061 : ICampbell fc Lattanziol 120081 ) . especially those that are still 
on the main sequence or sub-giant branch and whose 12 C/ 13 C ratios cannot therefore have been 
reduced by extra mixing during their own FGB. IStancliffd (|2010i ) has suggested that thermohaline 
mixing can explain the Li enrichment, but the models do not have low enough 12 C/ 13 C ratios. 

The oxygen isotope ratios measured in evolved giants also do not support a case for extra 
mixing on the AGB. The observed increase in the oxygen isotopic composition with evolutionary 
state can be explained by the action of the TDU and higher O intershell abundances than predicted 
by standard models (see also the discussion in iHerwid |2000| ) . Increasing O is motivated by both 
theoretical models of AGB stars as well as observations of hot post-AGB objects. One issue 
with this comparison are the large uncertainties in the measurements of oxygen isotope ratios in 
evolved stars. Updated measurements with smaller uncertainties would certainly help constrain this 
scenario. The oxygen isotope data for pre-solar oxide grains has been used as suppo r ting evidence 



for the existence of effic i ent extra mixing on the FGB and AGB (jZinner et al. 



2003 



Nittler et al 



1997 



2005 



Nollett et al 



Wasserburg et al.lll995l ). However, the most extreme grains fall outside of 



the region determined for stellar spectra. If these grains have an AGB origin, then either they are 
formed in stars too dusty and heavily obscured for spectroscopic analysis, or the stars that form 
such grains are rare and not included in observational samples. One test that would help to confirm 
(or rule out) such mixing would be to derive accurate measurements of oxygen isotope ratios for 
large samples of O-rich and C-rich AGB stars covering a range of metallicities. This is desperately 
needed if we are to sort out the conflicting results for oxygen, especially given the consistency found 
for carbon between pre-solar SiC grains and AGB stars. 

The models and tests we performed have many uncertainties, including the depth and efficiency 
of the third dredge-up, which may be stochastic in nature, as well as mass loss, opacities, and 
reaction rates. We have circumvented the uncertainties related to the origin of the extra mixing in 
FGB stars by adopting the observed abundance ratios at the tip of the FGB for our AGB models. 
We have also simply assumed that all stars experience extra mixing on the F GB such that the value 
of th e 12 C/ 13 C ratio after mixing is equal to 10, whereas the observations by lGilroy &: Brownl (e.g., 
199ll ) show a variation between ~ 10 to 25. Our results indicate that the observational data for C/O 
and 12 C/ 13 C could in fact be fit better by assuming values between these extremes. Furthermore, 
the question of AGB intershell abundances may be affected by some of these uncertainties. Is 
convective overshoot the only way to obtain large O intershell abundances? Instead, could variations 
in key reaction rates (e.g., 12 C(a, 7) 16 0) allow us to synthesize more 16 in the intershell? Clearly 
further investigation into the intershell abundances of AGB stars is warranted. 
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Fig. 1. — Histogram of observed 12 C/ 13 C ratios of AGB stars from lLambert et al.l (jl986l ). Data for 
the J-type carbon stars have been excluded. The average value (excluding the J-type stars) is 58, 
shown by the solid gray vertical line. With the J-stars, the average is 50. The standard deviations 
are 25 and 17, with and without J-stars, respectively. The top panel shows model predictions from 
the 1.25M0 model and the bottom panel for the 1.8-Mq case. In the top panel, the dashed vertical 
line shows the 12 C/ 13 C ratio at the tip of the AGB for the N ov = 3 model, and the dot-dashed 
vertical line shows the result when the 12 C/ 13 C ratio is reduced at the tip of the FGB. The dotted 
line shows the ratio for the N ov = 4 model with extra mixing on the FGB. In the bottom panel, 
the dashed vertical line shows 12 C/ 13 C ratio at the tip of the AGB for the Nay = 1 model with no 
extra mixing, whereas the results including extra mixing are shown by the vertical dot-dashed line. 
The dotted line shows the ratio for the N ov = 2 with extra mixing on the FGB. 
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Fig. 2. — T h e C/ ratio versus C/ C. Observational data of Galactic C-star s from 



Lambert et al.l (jl986l ) are shown as solid circl es, data for NGC 184 6 from iLebzelter et al.l ((2008) 
as solid squares, and data for NGC 1978 from iLederer et al.l (|2009h as solid stars. Error bars are 
included for the two Magellanic Cloud clusters, except for one star which has a lower limit of 200 
for 12 C/ 13 C. Open symbols (connected by various line types) indicate model predictions from the 
1.8M© model with N ov = 3 after each TDU episode. Predictions are shown for two models with 
a scaled solar composition, with and without extra mixing, and for two a-enhanced models, with 
and without extra mixing. Labels indicate the C/0 and 12 C/ 13 C ratios at the tip of the FGB, with 
the standard model indicated by the solid line. 
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Fig. 3. — The 16 Q/ 18 Q r atio versus 16 0/ 17 Q . The solid points s how o bservational data from 



Harris fc Lambertl (|l984l ), lHarris et all (|l985al ). lHarris et all (jl985bl ). and lHarris et all (| 19871 ) for 



G and K giants, barium stars, MS and S-type stars, a nd carbon stars. Error bars are not shown 
but uncertainties are large; see e.g.. lHarris et al.l (|1987l ) for details. The dotted lines show the solar 
values for the oxygen isotope ratios, which are the initial values used in the stellar models. The 
open diamonds show the surface composition after each TDU from the 1.8M© model with N ov = 3 
assuming a standard intershell composition, and open squares an 16 intershell abundance of 20% 
(by mass). The open circles show results from the 1.25M model with iV ov = 4 assuming a standard 
intershell composition, and open triangles assuming an 16 intershell abundance of 20%. 
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Fig. 4. — The C/0 ratio versus 12 C/ 13 C. Filled symbols are the same as in Fig. [2l Open symbols 
indicate model predictions from the 1.8M model with iV ov = 3 with extra mixing on the FGB. 
In this plot we show the result of varying the intershell 12 C and 16 abundances. The dotted line 
shows the result of using the composition of PG1159-035 which has 48% C and 17% O. The dashed 
and dot-dashed lines shows the results of increasing the 16 abundance from the standard ~ 0.6% 
in the intershell to 10% (dashed) and 15% (dot-dashed), while keeping 12 C constant at 24%. 
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Table 1: Initial composition, post-FDU, tip-FGB for the extra-mixing case, and tip of the AGB 
abundance results for C, N, and O. 



Position 


N a 

JV OV 


Mix? b C/O c 


12 


c/ 13 c 


12 C/ 14 N 


14 N /15 N 


16Q/17Q 


16Q/18Q 






1.25M , 


z 


= 0.01 scaled-solar model. 






Initial 


- 




0.500 




89.0 


3.60 


272.0 


2680 


500 


post-FDU 


- 





0.387 




26.6 


1.48 


798.1 


1801 


587 


tip-FGB 


- 


Yes 


0.335 




10.0 


1.00 


960.2 


1801 


587 


tip-AGB 


3.0 


No 


1.077 




75.6 


4.01 


884.0 


1240 


610 


tip-AGB 


3.0 


Yes 


1.015 




33.3 


3.13 


1005 


1229 


607 


tip-AGB 


4.0 


No 


2.213 




160 


8.38 


876.9 


1225 


615 


tip-AGB 


4.0 


Yes 


2.165 




73.8 


6.87 


1039 


1223 


615 








1.8M , Z -- 


= 0.01 scaled-solar 


model. 






Initial 


- 




0.500 




89.0 


3.60 


272.0 


2680 


500 


post-FDU 






0.326 




22.1 


0.99 


1269 


331.0 


678 


tip-FGB 




Yes 


0.320 




10.0 


0.91 


1295 


331.0 


678 


tip-AGB 


1 ' 


No 


1.244 




88.1 


3.87 


1309 


325.8 


694 


tip-AGB 


1.0 


Yes 


1.252 




40.9 


3.66 


1261 


328.2 


698 


tip-AGB 


2.0 


No 


1.960 




142 


6.20 


1292 


331.0 


703 


tip-AGB 


2.0 


Yes 


1.966 




65.8 


6.31 


1265 


330.4 


703 


tip-AGB 


3.0 


No 


3.318 




246 


10.7 


1217 


338.9 


717 


tip-AGB 


3.0 


Yes 


3.318 




120 


10.4 


1253 


338.2 


718 








1.8M , 


Z -- 


= 0.01 Or 


-enhanced 


model. 






Initial 






0.350 




89.0 


3.60 


272.0 


3803 


707 


post-FDU 






0.230 




22.1 


0.99 


1269 


344.0 


962 


tip-FGB 




Yes 


0.226 




10.0 


0.91 


1293 


344.0 


962 


tip-AGB 


3.0 


No 


2.378 




246 


10.7 


1207 


345.3 


1000 


tip-AGB 


3.0 


Yes 


2.380 




120 


10.4 


1238 


345.3 


1000 



"This is the overshoot parameter that governs the efficiency of the TDU. 
6 This indicates if extra-mixing is assumed to occur on the first giant branch. 
C A11 abundance ratios are by number. 
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Table 2: AGB abundance results for C and O for models with various intershell compositions. The 
1.25M model has N ov = 4 while the 1.8M model has N ov = 3. 



12/~< a 16p> a 


c/o 6 


12 C/ 13 C 


16Q/17Q 


16Q/18Q 


1.25M , Z = 


■ 0.01 scaled-solar model. 




0.26 c 0.0078 c 


2.18 


73.8 


1225 


615 


0.26 0.20 


1.07 


73.8 


249 


1248 


0.60 0.20 


2.27 


157 


249 


1248 


1.8M , Z = 


0.01 scaled-solar model. 




0.24 c 0.0058 c 


3.33 


120 


338 


718 


0A8 d o.n d 


2.49 


231 


869 


1843 


0.24 0.10 


1.76 


120 


642 


1364 


0.24 0.15 


1.40 


120 


805 


1706 


0.60 0.20 


2.78 


287 


966 


2049 



"Intershell abundance (in mass fractions). 

6 A11 abundance ratios are by number. 

c Standard intershell composition. 

d The composition of the post-AGB star PG1159-035. 
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